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ABSTRACT 

We have developed a stellar wind model for OB supergiants to investigate the effects of ac- 
cretion from a clumpy wind on the luminosity and variability properties of High Mass X-ray 
Binaries. Assuming that the clumps are confined by ram pressure of the ambient gas and ex- 
ploring different distributions for their mass and radii, we computed the expected X-ray light 
curves in the framework of the Bondi-Hoyle accretion theory, modified to take into account 
the presence of clumps. The resulting variability properties are found to depend not only on 
the assumed orbital parameters but also on the wind characteristics. We have then applied this 
model to reproduce the X-ray light curves of three representative High Mass X-ray Binaries: 
two persistent supergiant systems (Vela X- 1 and 4U 1 700-377) and the Supergiant Fast X-ray 
Transient IGR Jl 1215— 5952. The model can reproduce well the observed light curves, but 
requiring in all cases an overall mass loss from the supergiant about a factor 3 — 10 smaller 
than the values inferred from UV lines studies that assume a homogeneous wind. 

Keywords: X-ray: individuals: VelaX-l/4U 1900-40, 4U 1700-377, IGR J11215-5952. 
stars: supergiants. 



1 INTRODUCTION 



A new class of massive X-ray binaries has been recognized in the 
, last few years, mainly thanks to observations carried out with the 
■ INTEGRAL satellite. They are transient X-ray sources associated to 
O or B supergiant stars and characterized by sh ort outbursts. Thes e 
Supergia nt Fast X—ray Transients (SFXTs) iSguer a et all feOOSh : 
iNegueruela et alj (|2()0d)] are remarkably different from the classi- 
cal High Mass X-ray Binaries (HMXBs) with supergiant compan- 
ions, that are bright persistent sources, and also differ from the Be 
transients for their optical companions and shorter outbursts. 

The outbursts of SFXTs involve a high dynamic range, span- 
ning 3 to 5 orders of magnitudes, from a quiescent luminosity of 



lO^'^ erg s ^ up to the peak luminosity of 10"* 



10^^ erg s" 



The outbursts typically last a few days and are composed of many 
short flares with duration of a few hours. Besides these bright out- 
bursts, the SFXTs can display a fain ter flaring activity with lumi- 
nosity Lx = 10^^ - 10^" erg s'^ iSidoli et aljlioosh . Different 
mec hanisms have been proposed to explain the SFXTs properties 
('see lSidolil ( l200^ and references therein for a recent review). 

ISidoli etal] ^2007^. proposed that SFXTs outbursts are due to 
the presence of an equatorial wind component, denser and slower 
than the spherically symmetric wind from the supergiant, and possi- 
bly inclined with respect to the orbital plane of the system. The en- 
hanced accretion rate occurring when the neutron star crosses this 
wind component can explain SFXTs showing periodic outbursts. 



such as IGR Jl 1215— 5952 and also other SFXTs, assuming differ- 
ent geometries for the outflowing equatorial wind. 

Another possibility involves the gate d mechanisms due to 
transi t ions across the ce ntrifugal barrier iOrebenev & SunvaevI 
l2007h . lBozzoetalJ ( l2008l) showed that a centrifugal or a magnetic 
barrier can explain the SFXTs properties only if the supergiant 
wind is inhomogeneous, and the accreting neutron star has a strong 
magnetic field (> lO'^* G) and a long spin period (Papin 10"^ s). 

lin'tZan3 ( l2005h proposed that the SFXTs flares are produced 
by accretion of clumps of matter from the compani on wind. In the 
framework of the clumpy wind mode l pro posed bvlOskinova et all 
( l2007h . lWalter & Zurita Herai j2007h and lNegueruela et alj i20m 
proposed that what distinguishes the SFXTs from the persistent 
HMXBs with supergiant companions is their different orbital sep- 
aration: in persistently bright sources the compact object orbits the 
companion at a small distance (<2 stellar radii) where there is a 
high number density of clumps, while the transient emission in 
SFXTs is produced by accretion of much rarer clumps present at 
larger distances. On the other hand, the m onitoring with Swift of 
a sample of four SFXTs jSidoli et alJboOSh has demonstrated that 
these sources accrete matter also outside the bright outbursts, so any 
model should also account for this important observational prop- 
erty. 

The winds of O and B type stars a re driven by the momen- 
tum transfer of the radiation field and ILucv & Soloiiioi] j 19701) 
showed that the dominant mechanism is the line scattering. The 
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an alytical form u lation for line-driven winds has been developed 
by ICastor et all ( 1 19751) (CAK). In the CAK theory the mass lost 
by the star is smoothly accelerated by the momentum transferred 
from the stellar continuum radiation, and forms a stationary and 
homogeneous wind. However, both observational evidence and the- 
oretical considerations indicate that the stellar winds are variable 
and non homogeneous. Changes in the UV line profiles, reveal- 
ing wind variabilit y, have been obs erved on time scales shorter 
than a day (.Kudritzki & Pulsl I200Q1 . The X-ray variability ob- 
served in 4U 1700-377, Vela X-1 and ot her HMXBs can b e 
explained in terms of win d inhomogeneity IWhite et"al] ( Il983h : 
iKrevkenbohm et alj ( l2008h l and further indications for the pres- 
ence of clumps come from X— ray spectroscopy. For example, the 
X— ray spectrum of Vela X— 1 during the eclipse phase shows re- 
combination lines produced by a hot ionized gas and fluorescent 
K- shell lin es produced by cool and dense gas of near-neutral ions 
jSako et a l. 1999). These authors proposed that the coexistence of 
highly ionized and near neutral ions can be explained with an in- 
homogeneous wind, where cool, dense clumps are embedded in a 
low er density, highly ioniz ed medium. 

lLucv&Whitel ( ll980h suggested that the wind acceleration is 
subject to a strong instability since small perturbations in the veloc- 
ity or density distribution grow with time producing a variable and 
strongly structured wind. The first time-dependent hydrodynami- 
cal simulations of u nstable line-driven wind s were performed by 
lOwocki e t al* (T98i) and, more r ecently, by ' Runacres & Owockil 
l l2002h and Runa cres & Owockil r2005). All these simulations 
show that the line-driven instability produces a highly structured 
wind, with reverse and forward shocks that compress the gas into 
clumps. Moreover, the shock heating can generate a hot in ter-clump 
mediu m into which the colder clumps are immersed jCarlberj 
[l983). 

Based on these considerations, we study in this work the ex- 
pected variability and X— ray luminosity properties of neutron stars 
accreting from a clumpy wind. In the next two Sections we describe 
our clumpy wind model and the assumptions for the mass accretion. 
In Section|4]we study the dependence of the parameters introduced 
in the previous sections. In Sections |5] |6] and |7] we compare the 
X-ray light curves predicted with our model to the observations of 
three high mass X-ray binaries, showing that it is possible to repro- 
duce well the observed flaring behaviors. 



2 CLUMPY STELLAR WIND PROPERTIES 

In our model, where the dynamical problem has not been treated, 
we assume that a fraction of the stellar wind is in form of clumps 
with a power law mass distribution 



(1) 



in the mass range A/a - M\^. The rate of clumps produced by the 
supergiant is related to the total mass loss rate A/tot by: 



• /Aftot , -1 
A'ci — — — — clumps s , 



(2) 



where / = Afci /A/tot is the fraction of mass lost in clumps and 
< A/ > is the average clump mass, which can be computed from 
Equation ([Til. 

Clu mps are driven radia lly outward by absorption of UV spec- 
tral lines l lCastoretal.ll 19751) . From spectroscopic observations of 



O stars, iLepine & MoffatI j2008l) suggest that clumps follow on av- 
erage the same velocity law of a smooth stellar wind. We can then 
assume the following clump velocity profile without solving the 
dynamical problem: 



v,i(r)^v^ (^1 „ 0.9983:^) 



(3) 



where Voo is the terminal wind velocity, i?oB is the radius 
of the supergiant, 0.9983 is a parameter which ensures that 
f (J^ob) ~ 10 km s~^, and / 3 is a constant in the ra nge ^0.5-1.5 
ILamers & Cassinellil ( ll999l) : lKudritzki et alJ(ll989h l. 

Assuming that the clumps are confined by the ram pressure of 
the ambient gas, th eir size can be derived by the balance pressure 
equation. FoUowing fLucv & White ( 1980.) and .Howk et al, ( ,2000b , 
the average density of a clump is: 



pel = Pw(r) 



(4) 



where pvjij') is the density profile of the homogeneous (inter- 
clump) wind, ttw and are the inter-clump wind and the clump 
thermal velocity, respectively: ai, = Jsl^ and ai = J£2i_. 
is the Boltzmann constant, r„ and Tc are the temperatures of the 
inter-clump wind and of the clumps, respectively, and /j, is the mean 
atomic weight. The constant Cp = 0.29 accounts for the confining 
effect of the bow shock pro duced by the ram pressure around the 
clump ( lLucv&Whi3l980h . uj is the relative velocity between the 
wind and the clump (w = iiw — fci). Adopting ~ SxlO^cms"'^, 



Tc ~ lO'^ K, Tw ~ 10'' K and ^ = 1.3, we obtaiiB 



200 



(5) 



Since the density radial profile Pv,{r) of the homogeneous 
inter-clump wind (obtained from the continuity equation M — 

{r)v{r) — constant) is: 

r-gu(ro) 



Pw{r) = Pw(ro)- 



(6) 



where ro is a generic distance from the supergiant, from Equations 
@, IDl, we obtain: 



Pcl(r) = Pcl(ro) 



(7) 



where Pci(ro) = Pw('"o) x (a^ + Cf^co^)/ a^. lBouret et alj ( |200"^ 
analyzed the far-ultraviolet spectrum of 0-type supergiants and 
found that clumping starts deep in the wind, just above the sonic 
point Rs, at velocity v{Rs) ~ 30 km s^^. In the CAK model the 
sonic point is defined as the point where the wind velocity is equal 
to the sound speed (v{Rs) — Cs): 

0.9983i?oB 



R. 



(8) 



1 - 

Adopting typical parameters for O supergiants, from Equation ^ 
we obtain that the clumping phenomenon starts close to the photo- 
sphere (Rs ~ Rob)- Assuming spherical geometry for the clumps 
and that the mass of each clump is conserved, it is possible to obtain 
the expansion law of the clumps from Equation (O, with ro — Rs'. 

1/3 



Rciir) = Rci{Rs) 



'«ci(r) 



RMR= 



(9) 



^ This value is just an estimate of a typical value of ratio of clump to am- 
bient density, and will not adopted throughout. 
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From Equation l|9j we find that the clump size increases with the 
distance from the supergiant star. 

For the initial clump dimensions we tried two different distri- 
butions, a power law 



Nm^i oc clumps s 

and a truncated gaussian function: 

Nm^i oc ^=e ^ V " I clumps s 



(10) 



(11) 



i?ci,min)/(2A'^o-), and Na is a free parame- 



where cr = (i?ci,max 
ter. 

For any given mass clump we derived the minimum and max- 
imum values for the initial radii as follows. The minimum radius 
is that below which the clump is optically thick for the UV reso- 
nance lines. In this case gravity dominates over the radiative force 
causing the clump to fall back onto the supergiant. The momentum 
equation of a radiatively driven clump is: 



dr 



GMoE 



+ 9c + 9L 



(12) 



where Mob is the mass of the supergiant, is the radiative ac- 
celeration due to the continuum opacity by electron scattering, 
is the radiative acceleration due to line scattering. While Equa- 
tion l ll2b is an approximation that ignores the pressure gradient, 
the solution of the momentum equation differs only s lightly from 
the accurate solution derived bv lKudritzki et al.l l ll989l) : another as- 
sumption is tha t the photosphere has been treated as a point source 
jLamers & Cassi n ellL 199ft) . The radiative acceleration due to elec- 
tron scattering is: 



9c (, 



ac{r)LoB 



(13) 



where Ue is the opacity for electron scattering, and is given by (Jo — 
(jtj^, where ar is the Thomson cross section, is the number 
density of electrons, pci is the densi ty of the clump and - ^ob is 
the luminosity of the OB supergiant. iLamers & Cassinellil j 19991) 
found 0.28 < CTc < 0.35 cm^ g~^. Assuming a constant degree of 
ionization in the wind both o"e and go are constant. 

The radiative acceleration due to the line scattering is: 



kt~ 



(14) 



a'i?^ = 0.325 cm^ g~^, fe, a, 5 are the force multiplier pa- 
rameters, which are obtained with the calculation of the line radia- 



tive force (lAbbotllll982h. t is the dimensionless optic al depth pa- 
ramete r dCastor et al.l 19751) . According to the model of lHowk et al.l 
( |200Q I). we assume no velocity gradient inside the clump (of size 
I « 2_Rci), then we utilize the dimensionless optical depth parame- 
ter for a static atmosphere t = ae [ pdr. W(r) is the geometrical 
dilution factor jLamers & Cassinellil 19991) . given by: 



W{r) 



Rob 



(15) 



According to the Equation 0, we find that the number density of 
electrons in each clump is given by: 

RsviRs) 



ne(r) = no- „ , , 
where no = ne(-Rs) such that: 

Pcl(-Rs) 



no 



(16) 



(17) 



From Equation l ll2b we obtain that the minimum radius of the 
clump is given by: 

5L + (7c-3g = (18) 
where g^ is the acceleration due to gravity. Approximating t as 

(19) 



Mci 3 Mci 

Vci 2 TvR^ 



and assuming the force multiplier parameters calculated by 
IShimada etid] ( ll994 . with v{Rs) = 30 km s"^ and ro = Rs, 
we finally obtain the lower-limit for the clump radius: 

A - B\ l/(3<5-2a) 



-^cl , m 

where 



3AfcilO" 



B = 



C = 



AnfiemnWiRs 



(20) 

(21) 
(22) 
(23) 



We found t hat for the interesting range of the clump masses, the 
drag force jLucv & W hite 1980) values are less than 3% of the 
forces resulting from Equations l ll3t and l ll4t . Thus, in the deter- 
mination of the minimum clump radius, we can neglect this contri- 
bution. 

The upper-limit to the clump radius is obtained from the defi- 
nition of a clump as an over-density with respect to the inter-clump 
smooth wind: 



Pel ^ 



4ttRIv^{R, 



(24) 



where Mw is the mass loss rate of the homogeneous wind compo- 
nent (inter-clump). From Equation l l24t we obtain the upper-limit 
for the clump radius 



Rc 



^3MciE|wwCRs 



1/3 



(25) 



The lower and the upper limits for the clump radius (Equations 
1201 and [24l( depend on the supergiant parameter Tefj, Lqb, Rob, 
Voo, P, M^. Thus supergiants of different spectral type have, for 
any given mass of the clump, different minimum and maximum 
values for the initial radii distribution. However, these differences 
are smaller than a factor 10 for OB supergiants, and the intersection 
of the upper-limit and lower-limit functions ranges from Md ~ 



lO^^gto; 



10^ 



3 LUMINOSITY COMPUTATION 



The Bondi-Hoyle-Lyttl e ton a ccretion theory iHovle & Lvttleto j 
( Il939l) : iBondi & Hovld Il944}) 1 is usually applied to the HMXBs 
where a OB supergiant loses mass in the form of a fast stellar wind, 
(terminal velocity, Voo ~ 1000 — 2000 km s~^), that is assumed to 
be homogeneous and spherically symmetric. Only matter within a 
distance smaller than the accretion radius (Ra.) is accreted: 



i?a(r) 



2G'Mns 



(26) 
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Figure 1. Schematic representation of our clumpy wind model, d is the distance between the centre of the clump and the centre of accreting compact object. 
i?a is the accretion radius. 



where Mns is the neutron star mass and v^ciif) is the relative ve- 
locity between the neutron star and the wind. The fraction of stellar 
wind gravitationally captured by the neutron star is given by: 

Maccr = /5wind(r)Wrcl(r)7r7?a (27) 

where Pwindl*") is the density of the wind jPavidson & Ostriked 
1 1973b . 

The Bondi-Hoyle-Lyttleton accretion theory, that is based on 
a homogeneous wind, requires an important modification to prop- 
erly take into account the presence of inhomogeneity in a clumpy 
wind. In the homogeneous case the wind particles are deflected by 
the gravitational field of the neutron star, and collide with the par- 
ticles having the symmetric trajectory in a cylindrical region with 
axis along the relative wind direction. The collisions dissipate the 
kinetic energy perpendicular to this axis, and only the particles with 
a parallel kinetic energy component lower than the gravitational po- 
tential energy are accreted. The application of this accretion mech- 
anism to an inhomogeneus wind can lead to a partial accretion of 
the clump: when the distance d between the neutron star and the 
projection of the centre of the clump on the accretion cross section 
is smaller than the clump radius and d 7^ 0, only a fraction of the 
mass of the clump will be accreted (see Figure[T). In particular, if 
an incoming clump that crosses the accretion cross-section 7ri?a, 
is smaller than the accretion radius (i?ci < Ra) and d < Rd, the 
accretion cross-section in Equation in\ must be replaced by 



(28) 



and the mass accretion rate is A/accr ~ Pci{r)vrci{r) x E (see 
Figure[2lA). If ^ R^. we have three cases: 

• when Rci — d^ R^, (Figure|2lB), the accretion cross-section 
is given by: 

^2 



E = -KRi 



(29) 



• when Rd - d< R^ and R^ - d? < -Ra (FigureHC), the 
accretion cross-section is given by Equation i28t . 



• when - d < R^ and y/R^~^ > R^ (Figure|2]D), the 
accretion cross-section is given by: 

- x^dx — {xi. 



E = 4 



d)yi 



-Ra 



(30) 



The number density of clumps rid obeys the equation of con- 
tinuity TVci — 47vr^nci{r)vci{r), where iVci is the rate of clumps 
emitted by the OB supergiant (see Equation[TJ. Thus: 



nc\{r) 



iVcl 



clumps cm 



(31) 



47rr^i'ci(r) 

The rate of clumps accreted by the neutron star is given by: 

iVaccr = nci(r)uroi(r) X {nRi) clumps s~^ (32) 

When the neutron star accretes only the inter-clump wind, the 
X— ray luminosity variations are due to changes in its distance from 
the OB companion due to orbital eccentricity. The corresponding 
luminosity (see Equations 1261 and |27ll is: 



vind(<?!') 



Rns 



Mac 



Rm [{vL{r)+c^}^/^ 



(33) 



where (f) is the orbital phase, and Pv{r) is given by the Equation 
©. 

When the neutron star accretes a clump its X— ray luminosity 
is given by: 

r tA\ CMns Mci 

ix.cl(0 = -5 i „o ^rcl X E (34) 

Rns |nRl^ 

where Rci{r) is given by the Equation l[9]l, and E by the Equations 
l l28t - l l3Qt . If two or more clumps are accreted at the same time, the 
X— ray luminosity at the peak of the flare produced by the accretion 
is given by the sum of the luminosities (Equation |34t produced by 
the accretion of each clump. 
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Figure 2. Schematic view of different possibilities for the accretion of clumps of radius R^i. A: R^i ^ /?a and d (Eg. I28L B: R^i > and 
R^l-d^ i?a (Eq.[29j; C: R^i - d < R^ and ^JW~cP < R^ (Eq.|2U; D: R^i - d < R^ and y/R^^ - d'^ > R^ (Eq.[30). 




18 20 
log,o(Mci) (g) 



Figure 3. Upper (dashed line) and lower-limit (solid line) for the clump 
radius at r = Re. They have been obtained from Equations (20} and {25\ 
assuming the following parameters for the supergiant: Mqb = 30 Mq, 
Rob = 23.8R0,Mtot = 10"'^ M© yr-\Doo = 1700 km s-\ /3 = 1, 
vo = 10 km s-i, Afei/A/„ind = 0.7. 



4 APPLICATION OF THE CLUMPY WIND MODEL 

In this section we investigate how the X-ray luminosity and vari- 
ability properties depend on the different clumpy wind parameters 
and orbital configurations. 

As an example, we consider a binary con sisting of an 08.5 1 
star with Mqb = 30 M©, 7?ob = 23.8 R© jVacca et al.lll996h . 
and a neutron star with A/ns ~ 1-4 Mq, i?NS = 10 km. The 
parameters for t he supergiant w ind are the following: Mtot = 
10"^ Mq yr-i jPuls et al.|[T996h . = 1700 km s-\ /3 = 1, 
Vo — 10 km s~^, and the force multi plier parameters are k = 
0.375, Q = 0.522, and S = 0.099 (Shima da et aT]|l994l) . The 
corresponding upper and lower limits to the clump radius, derived 
from Equations J20| l and l |25l >. are shown in Figure [3] 



4.1 The effect of the mass distribution 

To study the effects of the clump masses we computed the distri- 
butions of the flares luminosity and durations for different values 
of f, A/a and Mh considering for simplicity a circular orbit with 
^orb = 10 d. We first neglected the clump radii distribution assum- 
ing that clumps which start from the sonic radius have radii given 
by Equation l l20t . and follow the expansion law (Equation |9}. We 



found that when A/a and/or Mh increase, the number of clumps 
produced by the supergiant decreases (see Equation^, resulting in 
a smaller number of X— ray flares. The average flare luminosity, the 
average flare duration, the number of flares and the shapes of the 
luminosity and flare duration distributions do not change much for 
different values of A/a and A/b. In Figure |4] we show the depen- 
dence of the distributions of flare luminosities and durations on 
and / (for A/a = 10^^ g, A/b = 10^^ g). If C increases, the number 
of clumps and their density increases (see Equations[2land|20l(. im- 
plying a higher number of flares, a shift to higher luminosities and 
to shorter flare durations. 

Figure [5] shows the effect of changing the fraction of wind 
mass in the form of clumps / for different values of C^. When / 
increases, the supergiant produces more clumps (see Equations |2](, 
thus the number of flares and their average luminosity increase, 
while their average duration remains unchanged. 



4.2 The effect of the radii distribution 

In Section|2]we showed that, for any given mass, the clump dimen- 
sions are constrained within the limits given by Equations ( 120b and 
( 1251 ). Here we show the effect of different assumptions on the radii 
distributions laws. We considered both a power law (Equation 1 10b 
and a truncated normal distribution (Equationlllb. described by the 
parameters 7 (or Ncr)- When 7 increases, the number of clumps 
with larger density decreases, thus the average luminosity of the 
flares decreases and their average duration increases (see Figure|6ll. 
We also found that when ( increases, the flare distributions with a 
positive 7 follow a different behaviour than those with a negative 
value: for 7 < 0, the flare distributions behave as described above 
(i.e. the number of flares increases with Q, while this does not hap- 
pen for 7 > 0. This is due to the fact that in this case the clumps 
are larger, thus there is a high probability that two or more clumps 
overlap thus reducing the number of flares. 

For the case of a normal distribution of clump radii, we calcu- 
lated the distributions of flare luminosities and durations for differ- 
ent values of Ncr (see Figure |7}. When Ncr increases, the number 
of clumps with larger and smaller radii is reduced, resulting in a 
narrower flare luminosity distribution. When ( increases, the flare 
distributions have the same behaviour described above for the case 
7 > 0. 

We then tried another test, increasing the orbital period (e.g. 
from 10 days to 100 days), finding that the shape of the integral 
distributions in Figure |6]remains similar, except for the number of 
X— ray flares, which decreases. 



6 L. Ducci et al. 



1500 



■-■ 1000 



CP 

o 



500 



34 



2500 
2000 
1500 
1000 
500 



34 



M„=10'^g, Mb=10^^g, f = 0.5 



M<,= 10'^g, Mb=10"g, f=0.5 



22, 



= 1.1 
= 1.5 
=2 
=2.5 



35 




36 37 
log(L) (erg/s) 



M„=10'^g, Mb=10^^g, f = 0.7 



35 



36 37 
log(L) (erg/s) 



38 



M„=10'^g, Mb=10^^g, f=0.9 



35 36 37 38 
log(L) (erg/s) 



39 





' ' ' ' 1 1 ' 


1 1 




f=l.l 






<-=1.5 


/ 




t=Z 


/ 

/ 




<-=2.5 


/ 
/ 






/ 
/ 

/ 

/ 






/ 

/ r 

i * 






ft ~ 






/ ■ 





















39 




39 



< 

V 



1500 



1000 



500 




10' 



1 

y* 

/ 

/ 

/ 


1 

^=2 

<-=2.5 


~ / 
/ 
/' 

i 




i 

i ' - ' ' 








~ ,' / 




" f / 




■ i / 

■ i / 

if' 





10 



10* 

At (s) 



2500 
2000 
1500 
1000 

500 




M<,= 10'^g, M(,= 102^g, f=0.7 



10' 



< 

V 



10= 



1 


1 


y 

i 


<-=1.5 

<-=2 - 

<-=2.5 : 


i 




- i 
i 
/ 




- ! 
■ / 
■ 

- 1 / 

:/ /' 




'.1/ '■^^^ 





10' 10* 
At (s) 

M„=10'^g, M(,= 10^^g, f=0.9 

-i-i-l r 

^=1.5 
<-=2.5 



10^ 
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4.3 The effect of the mass-loss rate 

The effects of different wind mass-loss rates Altot are shown 
in Figure [8] The mass-loss rate is usually derived observation- 
ally from the strength of the Ha emission line, since this gives 
smaller uncertainties th an the method based on UV P-Cygni lines 
jKudritzki & PulsllOOOl) . Since the Ha line opacity depends on p , 



the presence of wind inhomogeneities leads to an over-estimate 
of the mass-loss rate. In particular, the mass-loss rates from O 
stars derived from smooth- wind models measurements with the Ha 
met hod need to be reduced b y a f actor 3 to 10 if t he wind is clumpy 
[see Lepine & Moffat! ( |2008[) and lHamann et alj (^8)]. In Figure 
[8] we show that if Mtot decreases, also the number of flares de- 
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creases due to the reduced number of clumps. Also the average 
luminosity of the flares is reduced because the number density of 
clumps decreases resulting in a smaller probability that two or more 
clumps overlap. 



4.4 The effect of the orbital parameters 

In Figure |9] we show the effect of changing the orbital period Pa^h 
and the eccentricity e. We assumed — 1.1, 7 = —1.5, / — 0.7. 
When the orbital period increases, the number of flares emitted by 
the neutron star decreases (see Equation 131 1) and the neutron star 
accretes clumps with a smaller density (see Equation[9}, implying a 
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shift to lower luminosities and higher flare durations (see Figure[9l(. 
When the eccentricity increases, the neutron star accretes clumps 
with a higher density range (in general clumps are denser when 
they are closer to the supergiant), thus the luminosity range of the 
flares increases, as shown in Figure|9] 



5 COMPARISON WITH THE HMXB VELA X 1 

Vela X— 1 (4U 1900—40) is a bright eclipsing X-ray binary 
(Porb= 8.9 64 d, e ^ 0.09) formed by the B0.5 lb supergiant 
HP 77581 ([Brucato&KristiaiJll972l) (M = 23 M©, R= 30 Rq, 
dvan Kerkwiik et alj[T99l) ) and a pulsar with spin period ~ 283 s 
and mass 1.9 Mq. This source shows significant X-ray variability 
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Figure 8. Expected integral distributions of the flare luminosities and durations for different values of the mass-loss rate and for an assumed orbital period of 
10 days and a circular orbit, f = 1.1, 7 = —1.5, / = 0.7. The other binary system parameters are in the caption of Figure|4] The time interval coiresponds 
to 100 days. 



on short time-scales, with flares lasting from ~ 500 s to ~ 40000 s 
iHaberll (ll994V Krevkenbohm et al.' feOOS')]. 

Recentlv .Krevkenbohm et alj l i2008i') analyzed INTEGRAL ob- 
servations of Vela X— 1 obtained during a phase of high flaring ac- 
tivity, finding two kinds of flares: brief and bright flares softer than 
longer flares. They also found several off-stat es, during which the 
source is not detected (at least by /A^rEG^D. lKrevkenbohm et al.l 
1 I2OO8I) proposed that the short flares are caused by the flip-flop in- 
stability, while the long ones are due to the accretion of clumps 
ejected by the supergiant. The off-states are explained as due to the 
onset of the propeller effect when the neutron star crosses the lower 
density inter-clump medium. 

In this Section we apply our wi nd model to Vela X -1, as- 
suming the wind parameters derived bv lSearle etal] ( l2008h and re- 
ported in Table [T] 

In Figure [To] we compare the light curve measured with the 
ASM/RXTE instrument with that calculated with our clumpy wind 
model assuming a spherical symmetry for the outflowing wind. The 
ASMJRXTE count rate, measured in the 2 — 10 keV range, has been 
converted to the 1 — 10 keV luminosity using the average spectral 
param eters obtained by Orlandini et al.l ( Il99&) and the distance of 
2kpc dSadakane et al.ll985l) . The observed light curve of Vela X— 1 
is well reproduced by our clumpy wind model for Mtot = 2.1 x 
lO"'^ M0 yr^\ A'/a = 5 X 10^** g and A/b = 5 x 10^^ g, C = 1.1, 
/ = 0.75 and 7 = — 1. Acceptable light curves were also obtained 
with 7 = 1 and Mtot = 4 x lO"'^ Mq yr^^, and with a normal 
distribution law for the clump radii, with = 5. We point out that 
the average luminosity observed by ASM/RXTE out of the flares is 
~ 3 — 4 X 10''^ erg s~^. This luminosity is obtained in our model 
with the accretion of numerous clumps with low density. 



6 COMPARISON WITH THE HMXB 4U 1700-377 

4U 1700-377 l ljones et al.|[l973l) is a 3.412 day eclipsing HMXB 
composed of a compact object, (a neutron star or a black hole), 
and the 06.5 laf'' s tar HD 153919, located at a distance of 1.9 kpc 
jAnkav et al]|200lh . Despite extensive searches, no X— ray pulsa- 
tions have been found in this system. Therefore the X-ray mass 



function cannot be determined and the system parameters {Mqb, 
Rob, A/a:,) cannot be derived directly. They have been estimated 
from the radial velocity curve of the supergiant and from the dura- 
tion of the X-ray eclipse, by making assumptions about possible 
values of the radius of the O star and the orbital inclination. Several 
studies indicate that the m a ss of the compact object is larger than 
2 Mq iRubin et alj ( Il996h : IClark et alj ilOOj) ]. The similarity of 
the X— ray spectrum to other pulsars s uggest that the com pact ob- 
ject of 4U 1700-377 is a neutron star jwhite et al.ll 1983b, but the 
prese n ce of a low-ma ss black hole cannot be excluded ( Brown et"al] 
fl993) .l Revnolds et al. ( 1999) reported the presence of a possible 
cyclotron feature at ~ 37 keV observed with BeppoSAX. If con- 
firmed, this would demonstrate that 4U 1700—377 hosts a neutron 
star with a magnetic field of about 2.3 x 10^^ G. 

The X— ray light curve of 4U 1700—377 is characterized 
by a strong flaring activity with variations as larg e as a factor 
of 10 — 100 on time scal e s from minu t es to hours iHaberl et al.l 
l l 19891) ; I White et al.l ( ll983h : lRubin et all (Il996^1. We assumed the 
most recent set of system parameters of 4U 1700—377, obtained 
by means of a detailed NLTE (Non-Local Thermal Equilibrium) 
line-driven wind model analysis of HD 153919 and a Monte Carlo 
simulation for th e determination of the masses of both components 
dClark et al. 2002). These authors found that the supergiant has a 
luminosity log(L/LQ) = 5.82 ± 0.07, an effective temperature 
Tcft ^ 35000 K, radius Rqb « 21.9 R©, mass Mqb ^ 58 Mq, 
mass loss rate M = 9.5 x 10~® Mq yr~^, and a mass for the 
compact object A'/x = 2.44 Mq (Table[T]l. 

The X-ray spectrum of 4U 1700—377 is well described by an 
absorbed power law with high-energy cutoff dvan der Meer et al] 
I2OO5I) . The spectrum above 20 keV was studied using different 
satellites (e.g. BATSE detector on board the CGRO, INTEGRAL) 
and can be mo delled using a ther mal bremsstrahlung model with 
kT ~ 25 k eV jRubin et alj|l996h or with a thermal Comptoniza- 
tion model jOrr et al .II2OO4I) . 

The analysis of the 0.5 — 12 keV spectrum with XMM- 
Newton during the eclipse , the egress, and a low-flux interval led 
Ivan der Meer et a i] ( l2005h to suggest that the low-flux interval is 
probably due to a lack of accretion such as expected in a str u ctured 
and inhomogeneous wind. Moreover, Ivan der Meer et all ( l2005h 
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Figure 9. Expected integral distributions of the flare luminosities for different orbital periods and eccentrities of the binary system. ^ = 1.1, 7 = —1.5 and 
/ = 0.7. The other binary system parameters are in the caption of Figure|4] The time interval for each histogram corresponds to the orbital period. 



proposed that the fluorescence line from near-neutral iron detected 
in all spectra is produced by dense clumps. They also observed re- 
combination lines during the eclipse which indicate the presence of 
ionized zone around the compact object. 

We analyzed the public archival INTEGRAL data of 
4U 1700-377, using all the IBIS/ISGRI observations obtained 



from 2003 March 12 to 2003 April 22, and from 2004 February 
2 to 2004 March 1. These data correspond to a net exposure time 
of ~ 5.2 days (excluding the eclipse phase). We reduced the data 
using OSA 7.0, and extracted the light curve in the energy range 
15—60 keV, finding a total of 123 flares. For each flare we extracted 
the spectrum in the range 22 — 100 keV. All the spectra could be 
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Table 1. Parameters of the HMXBs studied in tliis paper. Values labeled with (*) are taken from lSearle et alj i2008l) , values labeled with (t) are taken from 
iLefever et all )2007l) . 





VelaX-1 


4U 1700-377 


IGRJ11215-5952 


Type 


SGXB 


SGXB 


SFXT 


Spectral Type 


B0.5 lb 


06.5 Iaf+ 


B0.7 la 


Mob 


23 Mq 


58 Mq 


23 Mq (*) 


Rob 


30 Rq 


21.9 Rq 


35.1 Rq (*) 




2.5 X 10-* K(*) 


3.5 X 10* K 


2.36 X 10* K(*) 


log(L/L0) 


5.58 (*) 


5.82 


5.5 


Mtot 


0.7 - 1.2 X 10-"^ Mq yr-l (*) 


9.5 X lO-'^ Mq yr-l 


1 - 2.5 X 10"^ Mq yr-l 




~ 1520 - 1600kms-i (*) 


1700kms-i 


~ 1000 - 1400 km s"! (*) (t) 




~ 1 - 1.5 (*) 


1.3 


~ 1 - 1.5 


Mns 


1.9 Mq 


2.44 Mq 


1.4Mq 


^orb 


8.964 d 


3.412 d 


164.5 d 
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e ~ 0.09 







P . 
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Figure 11. Comparison between observed (solid line) and calculated (dashed line) distributions of the flare luminosities and durations for 4U 1700—377. 
The binary system parameters are: A/qb = 58 Mq, Rqb = 21.9 Rq, M^s = 2.44 Mq, i?iMS = 10 km. The parameters for the supergiant wind are: 
Mtot = 1.3 X 10-^ Mq yr-\ Voo = 1700 km s-\ (3 = 1.3, vq = 10 km s"!. A/a = 5 X lO^^^ g and Aft = 2 X lO^^ g, C = 1-2, 7 = -6.5 and 
/ = 0.75. 



well fit by a thermal Comptonization model (COMPTT in XSPEC). 
Based on the spectral results, we computed the 1 — 200 keV lu- 
minosity of each flare. All of them have a luminosity greater than 
5.8 X 10"^^ erg s~^ (for lower luminosities it is difficult to evaluate 
the flare duration and then to distinguish the flares from the aver- 
age level of the X— ray emission). For each flare we have measured 
two parameters: the peak luminosity and the flare duration. We then 
applied our clumpy wind model to the INTEGRAL observations of 
4U 1700—377. We first compared the observed distributions of the 
flare luminosities and durations with those computed adopting the 
system parameters reported in Table [T] We assumed for the com- 
puted distribution a time interval equal to the exposure time of the 
4U 1700—377 observations considered here. As shown in Figure 
[TT]the flare properties are well reproduced with our clumpy wind 
model for A/tot = 1.3 x lO"** Mq yr-\ A/a = 5 x lO^** g and 
Mb = 2 X lO^'' g, C = 1-2, 7 = -6.5 and / = 0.75. We found 
that the numbers of observed (123) and calculated flares (116) are 
in good agreement. The light curve of 4U 1700—377 calculated 
with our clumpy wind model is shown in Figure[T2l 



7 COMPARISON WITH THE SFXT ICR J11215 5952 

The SFXT IGR J 112 15-5952 was di scovered in April 2005 
with INTEGRAL fcubinski et J] |2005|) . It is associated with 
HD 3 06414, a B0.7 la star lo cated at an estimated distance of ~ 

8 kpc jNegueruela et al.|[2007h . R XTE observa t ions s h owed a pulse 
period P^pin = 186.78 ± 0.3 s ISmith et alj ( l2006h : ISwank et al.1 
(12007)]. This is the first SFXT for whi ch a periodicity in the out- 
bursts recurrence time was discovered jSidoli et al.|[200^ . Subse- 
quent observations jRomano et alj|200^ showed that the true pe- 
riodicity is about 164.5 days, i.e. half of the originally proposed 
value. This periodicity is very likely due to the orbital period of the 
system. 

For a distance of 8 kpc the peak fluxes of the outbursts cor- 
respond to a luminosity of ~ 5 x W^^' erg s^^ [5 — 100 keV, 
Romano et al. ( 2007)] • Swift monitoring of this source revealed that 
the outburst (lasting a few days) is composed by many flares (last- 
ing from minutes to a few hours), and before and after the whol e 
outburst the source is fainter than 10^^ erg s~^ jSidoli et"al]|2007h . 

ISidoh et"ai] ilOOt showed that accretion from a spherically 
symmetric homogeneous wind could not reproduce the observed 
light curve and therefore proposed a model based on an anisotropic 
wind characterized by a denser and slower equatorial component. 
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Figure 10. Comparison of the Vela X— 1 light curve, as observed with 
ASMJRXTE (lower panel), with that calculated with our clumpy wind 
model for the parameters reported in Table [T] and Aftot = 2.1 X 
10-'^ Mq yr-l, (3 = l,Voo = 1600 km s"!, / = Md/M^^ind = 0.75. 
Ma = 5 X lOl* g, Mb = 5 X 10^1 g, C = 1.1, = 5 (middle panel), 
and 7 = -1 (upper panel), k = 0.375, a = 0.522, & = 0.099. Orbital 
phase = con'esponds to 53750 MJD. Note that Vela X— 1 is not contin- 
uously observed with ASMJRXTE. Therefore it is possible that some flares 
have been missed. 



37.4 



periodically crossed by the neutron star along its orbit. However 
this result was based on the old determination of the orbital period 
(329 days). Therefore, before applying our clumpy wind model, 
we checked the spherically symmetric homogeneous wind with 
Poib = 164.5 d, different eccentricity values, and t he set of stel- 
lar parameters derived by ISearle etal] | |2008|) and iLefever et al.l 
(2007) and reported in Table [T] We assumed a terminal velocity 
Voo ranging from 1000 km s~^ to 1400 km s~^, and Af ranging 
from 1 X 10"^ Mq yr"^ to 2.5 x 10"" Mq yr"\ In all cases we 
found that the duration of the X— ray outburst observed with Swift 
is shorter than that of the calculated light curves. 

A better agreement with the observations could be obtained 
with our clumpy wind model, especially for what concerns the flar- 
ing variability during the outburst phase. However, also in this case 
the calculated light curve always produces an outburst longer than 
the observed one. This is shown in Fig. [T3] where the green sym- 
bols corresponds to the accretion of a dense clump (producing a 
flare), while the blue symbols indicate the lower luminosity level 
produced by the accretion of the inter-clump matter. We found that 
the probability to observe a flare, rather than the inter-clump lumi- 
nosity level, is 90%. 

In order to improve the agreement between the observed and 
the calculated light curve, we introduce d an a nisotropic outflow 
similar to that proposed by ISidoli et al.l ilOOt . In this modified 
model we introduce a denser clumpy wind component in the equa- 
torial plane, with a thickness 2h, a terminal velocity v^d and a mass 
loss rate Al^d, together with a spherically symmetric clumpy wind 
component (polar wind) with terminal velocity Vpw and a mass loss 
rate A/p„ . We linked the mass loss rate from the equatorial outflow 
with the mass loss rate from the polar wind, by means of the factor 

fad- 



Me, 



/odMp 



(35) 



Both wind components are clumpy and obey laws described in Sec- 
tion |2] We assume that the second wind component has a Gaus- 
sian density profile perpendicular to the equatorial plane of the su- 
pergiant. In this framework we have considered an orbital period 
forb = 164.5 d and a high eccentricity in order to produce only one 
outburst per orbit. The comparison between the Swift light curve 
and that predicted with this model is shown in Figure [14] A good 
agreement with the data is obtained with the parameters reported in 
Tabled 




Figure 12. Computed hght curve of 4U 1700-377. The off-state at orbital 
phase 0.42 ^ 4> ^ 0.58 is due to the eclipse. The stellar parameters are 
reported in the caption of Figure [TT] 



8 DISCUSSION AND CONCLUSIONS 

We have developed a clumpy wind model (where the dynamical 
problem is not treated) and explored the resulting effects on an 
accreting compact object in order to explain the observed behav- 
ior of the SFXTs and the SGXBs. Compared to previous attempts 
to explain the SFXTs outbursts in the co ntext of clu mpy winds 
[Walter &ZuritaHeras (2007); Negueruela et alj ((2008)1. we in- 
troduced a distribution for the masses and initial dimensions of 
the clumps. We described the subsequent expansion of the clumps 
(Equation|9) taking into account realistic upper and lower limits for 
their radius (Eauationsl20landl25t. 

This model, together with the theory of wind accretion modi- 
fied because of the presence of clumps, allow us a comparison with 
the observed properties of both the light curves and luminosity dis- 
tributions of the flares in SGXBs and SFXTs. 

From the calculated integral distributions (Section |4ll, we 
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Figure 13. Comparison of the IGR Jl 1215— 5952 light curve of the February 2007 outburst observed with Swift/XKt with the calculated X— ray light curve of 
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Figure 14. Comparison of the IGR Jl 1215—5952 light curve of the February 2007 outburst observed with Swift/XKt with the calculated X— ray light curve 
of our clumpy wind model with equatorially enhanced wind component (presented in Section|2). We used the stellar parameters reported in Table|2] 



found that the observable characteristics of the flares, such as lu- 
minosity, duration, number of flares produced, depend mainly on 
the orbital period (Figure[9}, the scaling parameter ( of the power- 
law distribution for the clump formation rate (Equation [TJ, and the 
fraction of wind mass in the form of clumps (/ = Md/Mtot), as 



shown in Figures|4]and[5l Thus the variability properties of the dif- 
ferent systems do not depend only on the orbital parameters, but 
are also significantly affected by the properties of the clumps (in 
particular by the parameters 7, /, Ma, Mb). 

We successfully applied our clumpy wind model to three dif- 
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Table 2. Parameters for the calculated X— ray light curves of 
IGR Jl 1215— 5952 of our clumpy wind model with equatorially enhanced 
wind component (see Figure [T4l . io is the angle of inclination of the orbital 
plane with respect the equatorial wind component; ip is the angle between 
the orbital plane intersection with the equatorial wind component and the 
direction of the periastron. 



Parameter 



Value 



13 

Mpw 

"^oo, pw 

io 
ip 

f = M,i/M„ind 

Ma 

Mb 

c 

7 

fed 

^oo , ed 



1 

1 X 10~'^ Mq yr" 
1500kms-i 
20° 
70° 
2. X 10" cm 
0.75 
5 X 10^** g 
1020 g 
1.1 
6 

100 
lOOOkms-i 



for the 3 sources. This table shows the typical differences in the 
mass-loss rate and terminal velocity expected from the CAK the- 
ory, similar values for ( and /, and different values for 7, i\fa, A/b, 
which seem to be proportional to the inverse of the effective tem- 
perature of the supergiant. Moreover, we found that the values of 
A/tot that best reproduce the observed light curves of the 3 HMXBs 
studied, are in agreement with the hypothesis that the mass-loss rate 
derived by the Ha emission is overestimated by a factor 2-10, in 
agreement with recent studies (see Section \43[ . We can exclude 
that the observed light curves of the three HMXBs studied can be 
reproduced with the same set of wind parameters. 

In conclusion, the different values of 7, A/a, Mh obtained for 
the 3 HMXBs studied in this paper reveal that, in the framework 
of our clumpy wind model, the properties of the clumps of these 
3 X-ray binary systems are slightly different, independently of the 
orbital period. This discrepancy could be due to the different spec- 
tral type of the 3 supergiants, which could eject structurally inho- 
mogeneus winds with slightly different properties. We suggest as a 
possible cause of this behaviour that the values of 7, Ma, Mb could 
be rel ated to the sonic radius where the clumps start jBouret et al.l 
l2005l) . which depends on the supergiant properties (Equation[8j. 




-0.3 -0.2 -0.1 0.1 0.2 0.3 

Orbital Phase 



Figure 15. Calculated light curve of a SFXT with a neutron star accreting 
an anisotropic clumpy wind. The parameters of the system are: Mjms = 
1.4 Mq, Mob = 23 Mg, Rps = 35.1 R©, Lqb = 3.16 x 10^ Lq, 
Porb = 80 d, e = 0.75, Mtot = 2 x lO"'^ Mq yr-\ Vaa,pw = 
1500 km s-\ /3 = 1, io = 15°, ip = 0°, h = 2. x lO" cm, 
Voo,ed = 1000 km s-\ /ed = 50, Teff = 20000 K, fc = 0.709, 
a = 0.470,(5 = 0.089, C = 1.1,7 = 6, Ma = 5xl0i*g, Mb = lO^" g. 

ferent high mass X-ray binaries: Vela X— 1, 4U 1700—377 and 
IGR J11215— 5952. For the latter source, however, we had to intro- 
duce a denser equatorial component (still with a clumpy structure) 
in order to reproduce the flare duration. 

In Figure [TS] we show an example of a light curve of 
a generic SFXT calculated in the case of anisotropic clumpy 
wind. We assume for the generic SFXT properties similar to 
IGR J1I2I5— 5952, with an orbital period of 80 d and an eccen- 
tricity e = 0.75. The orbital plane intersects the equatorial wind 
component at two phases (cfii ~ —0.12, (j>2 ~ 0.12) producing two 
outbursts. The third outburst is produced at the periastron passage 
((f>3 ~ 0). This implies that, if this explanation is correct, up to 3 
outbursts per orbit are possible. 

In Table[3]we have summarized the wind parameters obtained 
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